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ABSTRACT 

Utilizing Chandra X-ray observations in the All-wavelength Extended Groth Strip International 
Survey (AEGIS) we identify 241 X-ray selected Active Galactic Nuclei (AGNs, L 2 -iokcV > 10 42 
ergs s _1 ) and study the properties of their host galaxies in the range 0.4 < z < 1.4. By making 
use of infrared photometry from Palomar Observatory and BRI imaging from the Canada-France- 
Hawaii Telescope, we estimate AGN host galaxy stellar masses and show that both stellar mass 
and photometric redshift estimates (where necessary) are robust to the possible contamination from 
AGNs in our X-ray selected sample. Accounting for the photometric and X-ray sensitivity limits of the 
survey, we construct the stellar mass function of X-ray selected AGN host galaxies and find that their 
abundance decreases by a factor of ~2 since z ~ 1, but remains roughly flat as a function of stellar 
mass. We compare the abundance of AGN hosts to the rate of star formation quenching observed in the 
total galaxy population. If the timescale for X-ray detectable AGN activity is roughly 0.5-1 Gyr — as 
suggested by black hole demographics and recent simulations — then we deduce that the inferred AGN 
"trigger" rate matches the star formation quenching rate, suggesting a link between these phenomena. 
However, given the large range of nuclear accretion rates we infer for the most massive and red hosts, 
X-ray selected AGNs may not be directly responsible for quenching star formation. 
Subject headings: cosmology: observations, galaxies: formation, galaxies: evolution 



1. INTRODUCTION 

Recent observations of the galaxy population and its 
evolution since z w 2 reveal a pattern in which the most 
massive galaxies appear to shut down star formation ac- 
tivity at early times with increasingly less massive galax- 
ies fo llowing later (e.g. , Uuneau et alJ 120051: rTreu et alJ 
20051; IBundv et alj2006tlBorch et al.ll2006UCimatti et all 



20061) . This pattern of "quenching" in the star forma- 
tion history of galaxies — thought to be largely respon- 
sible for the g rowing abundance of galaxies on the red 
sequence (e.g., iFaber et al.ll2007l) — is c ommonly referred 
to as "downsizing" (|Cowie et al.|ll996l) . 

Given that the dark matter halos of galaxies are 
expected to assemble hierarchically in the ACDM 
paradigm, understanding the physical mechanisms re- 
sponsible for downsizing remains an important challenge. 
We seek a process capable of quenching star forma- 
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tion, driving galaxies onto the red sequence, and pre- 
venting further star form ation. Work with the DEEP2 
Galaxy Redshift Survey (|Davis et alJ l2003f ) has shown 
that such a process must operate over a range of en- 
vironmental densities (IBundv et ail l2006t ICooper et al.1 
l2007t iGerke et al.ll2007l ). suggesting an internal compo- 
nent to quenching that acts in addition to the suppres- 
sion of star formation expected in high-density environ- 
ments. Meanwhile, theoretical work has focused on one 
such internal pr ocess, namely the p otential role played by 
AGN feedback ()Silk fc ReeJll998f ). both as a way of ex- 
plosively initiating the quenching event (iGranato et al.1 
12004 IScannapieco et aT]|2005t iHopkins et al.ll2005a|) and 
preventing hot ga s in already pa ssive galaxies from cool- 
ing to form stars (ICroton et al.ll200l lBower et al.l l2006l: 



ICattaneo et al.ll2006t Ide Lucia fc Blaizotj|2007i r 

These scenarios remain largely untested because ob- 
servational evidence has been difficult to obtain. One 
of the most promising ways forward is to examine the 
properties of AGN host galaxies and search for sig- 
natures of this feedback. Such observations are chal- 
lenging however because no sele ction method find s all 
galaxies that host active nuclei (|Mushotzkvl I2004T ) and 
brighter AGNs (quasars) can easily outshine their hosts. 
With such difficulties in mind, previous studies suggest 
that most lower luminosity AGNs t end to be found in 
massive, mostly spheroidal galaxies (Dunlop et al. 200 



Kauffmann et al 



spheroidal galaxies I Dunlo p et al.llz UUcS; 
1 l200l iGrogin et all 120051 : iPierce et all 



2007f ) although some "tran sition" sources exhib i t dis 



turbed morphologies fe.g.. ICanalizo fc Stocktonl 120011 : 
iHutchings et al.l 120061 : IConselice et al.l 120071 ). AGNs de- 



tected through emission line diagnostics and X-rays ap- 
pear to prefer host galaxies on th e red sequence and 
"green valley" (jNandra et all l2007t iMartin et all 120071 : 
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iSalim et al.ll200l . 

Further progress on testing the link between AGN ac- 
tivity and the downsizing of star formation requires that 
we understand the stellar mass distribution and redshift 
evolution of AGN hosts as compared to evolutionary pat- 
terns in the general galaxy population. In this paper we 
use deep Chandr a observations to i dentify AGN hosts in 
the AEGIS field (|Davis et al.ll2007l ) where spectroscopic 
and photometric redshifts (photo- z 's) as well as infrared 
photometry provide reliable stellar mass estimates out to 
z = 1.4. We then compute the AGN host stellar mass 
function and use it to estimate the rate at which AGN 
activity is triggered in the galaxy population. We will 
show that this mass-dependent rate is consistent with 
the rate of star formation quenching of all galaxies at 
the same epochs at which the AGN activity is observed. 

The structure of the paper is as follows. We describe 
the multi-wavelength observations and properties of the 
sample in $2] and Ej3j The way in which stellar masses 
and restframe colors are determined is given in §4] while 
our methods for constructing mass functions and the re- 
sulting AGN host mass function are presented in Sj5j In 
§6] we present evidence for a link between AGN activ- 
ity and quenching in the context of estimates of the X- 
ray AGN timescale and explore whether feedback from 
X-ray selected AGNs causes quenching. We summarize 
in |J71 Where necessary, we assume a standard cosmo- 
logical model with fijyi = 0.3, Qa = 0.7, H — 70/i7o 
km s" 1 Mpc" 1 . 

2. OBSERVATIONS 

In this section, we discuss the various observations we 
use to investigate the link between AGNs and galaxy evo- 
lution. We begi n by summarizing the DEEP2 Galaxy 
Redshift Survey (|Davis et alj I2003T ) and a followup in- 
frared imaging campaign conducted at Palomar Obser- 
vatory. One of the four fields in this survey, the Ex- 
tended Groth Strip (EGS), was imaged with Chandra as 
part of the AEGIS program, providing a sample of AGN 
host galaxies that we utilize here. We will compare the 
evolution of this AGN host sample to trends observed 
in galaxies dra wn from the f ull DEEP2/Palomar survey. 
As discussed in I Bund v et al.1 (|2006h . this sample provides 
robust stellar mass estimates that can be used to inves- 
tigate mass-dependent evolution out to z = 1.4. 

2.1. The DEEP2/Palomar Survey 

We provide a brief summary of the DEEP2 Galaxy 
Redshift survey and the near-infrared (IR) followup 
imaging conducted at Palom ar Observatory. Th e data 
sets involved were utilized by iBundv et al.1 (|2006l ) to ex- 
plore the nature of star formation downsizing since z « 
1.2, and further details are provided there. Now com- 
plete , the DEEP2 Galaxy Redshift survey (|Davis et al.l 
120031 ) utilized DEIMOS on the Keck-II telescope to ob- 
tain spectroscopic redshifts for ^40,000 galaxies with 
z < 1.5. The survey is magnitude limited at Rab < 24.1 
and covers more than 3 square degrees over four fields, 
one of which is the AEGIS field and is described further 
below. 

Redshift targets were selected using BRI photometry 
from the Canada-France-Ha waii Telescope (CFHT) with 
the 12Kx8K mosaic camera jCuillandre et al.ll2001h . As 
described in ICoil et all (|2004f ). the photometry reaches 



Rab ~ 25.5 and was also used for estimating photo- 
z's and restframe (U — B) colors, as discussed be- 
low. In the three non- AEGIS fields, target selection 
in DEEP2 was carried out using observed colors to ex- 
clude sources with z < 0.7. These selection criteria 
successfully recover 97% of the Rab < 24.1 popula- 
tion at z > 0.75 with o nly ~10% contam ination from 
lower redshift galaxies (| Davis etafl l200l . More de- 
tails on the observing strategy and character is tics o f 
the DEEP2 s a mple ar e provided i n ICoil et al.l (12004ft. 
Will mer et alj (120061). IDavis et alj (|2005h . iFaber et all 
(|2007f) . IDavis et al.l (|2007l ). and Newman et al. (in prepa- 
ration). 

Followup imaging of the DEEP 2 survey in the J and 
-K^-band were carried out using Wide Field Infrared 
Camera (WIRC. IWilson et all 2 00a ) on the 5m Hale Tele- 
scope at Palomar O bservatory. These o bservations are 
discussed in detail in lBundv et al.1 (|2006( ). Excluding the 
AEGIS field, which was the highest priority, 0.9 square 
degrees of DEEP 2 were imaged primarily in the JT s -band 
with exposure times varying from 2 to 8 hours, depending 
on the conditions, and the final typical 80% completeness 
depth of Kab = 21.5. 

We used SExtractor (jBertin k. Arnoutd 119961 ) to de- 
tect and measure 7^-band sources and cross-referenced 
them with the CFHT optical and DEEP2 redshift cat- 
alogs to construct a K-selected sample that forms the 
basis of our analysis. For total magnitudes used to esti- 
mate stellar masses, we took the MAG_AUTO output 
from SExtractor and did not correct this Kron-like mag- 
nitude for missing light. We estimated the uncertainty 
on these magnitudes by inserting fake sources at vari- 
ous magnitudes and using SExtractor to recover them. 
Based on the locations of K-selected sources, we mea- 
sured aperture photometry in the BRI data using 2" 
diameter apertures which we found exhibited the least 
scatter. These colors are used in fitting galaxy spectral 
energy distributions (SEDs) needed for both photo-z's 
and stellar mass estimates. 

2.2. The AEGIS Field 

Covering the 0.5 deg 2 Extended Groth Strip (a = 
14 h 17 m ,<5 = +52°30'), the All- Wavelengt h Extended 
Grot h Strip International Survey (AEGIS, IDavis et al.1 
|2007| ) accounts for one of the four fields where DEEP2 
redshifts and Palomar near-IR imaging were obtained, 
although there are slight differences with respect to the 
other DEEP2 fields (above) that we describe here. This 
field is also special because many other observatories, 
including the Hubble and Spitzer Space Telescopes, have 
conducted observations there. Of key importance for this 
work are the deep X-ray data from the Advanced CCD 
Imaging Spectrometer (ACIS) on Ch andra. 

The Chan dra data are describ e d inlGeorgakakis et al.1 
(2006) and iGeorgakakis et all fl 2007b]), and t he X - 
ray data analysis is presented in iNandra et al.1 (2005). 
Briefly, the AEGIS region was targeted over several 
epochs for a total integration time of 190 ks. Standard 
reduction methods using the CIAO software were em- 
ployed to derive fluxes in four energy bands — 0.5-7.0 
keV (full), 0.5-2.0 keV (soft), 2.0-7.0 keV (hard), and 
4.0-7.0 (ultrahard) — by integrating the counts detected 
within the 70% encircled energy radius at each source po- 
sition. The counts in each observed band were converted 
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to the standard bands of 0.5-10, 0.5-2, 2-10, and 5-10 
keV by assuming an intrinsic power-law with T = 1.4 
and Galactic absorption. The typical detection limits in 
each band are 35, 1.1, 8.2, and 14 in units of 10~ 16 ergs 
s _1 cm" 2 , although these depend on position because the 
Chandra sensitivity declines away from the center of each 
pointing. X-ray detections were referenced to the opti- 
ca l/IR photometric catalog s using the method described 
in IGeorgakakis et alj (|2006l ). 

In terms of the DEEP2 spectroscopic redshifts, target 
selection was carried out differently in AEGIS than it 
was in the other DEEP2 fields. Here, a straight mag- 
nitude limit of Rab < 24.1 was employed with galaxies 
having z < 0.7 down- weighted but not excluded from the 
target sample. The near-IR Palomar imaging is deepest 
in EGS, especially along the center of the strip where 
J-band imaging was also obtained, and the typical 80% 
completeness depth is Kab — 22.5. Overall, the typical 
AEGIS depth reaches Kab ~ 22. The properties of the 
CFHT BRI data are identical to those described above. 

3. THE GALAXY AND AGN HOST SAMPLES 

To begin our analysis of AGN hosts and their rela- 
tionship to evolution in the galaxy population we dis- 
tinguish between the galaxy sample and the AGN host 
sample. The galaxy sample is comprised of the full 
DEEP2/Palomar survey, including the AEGIS. With the 
goal of providing a self-consistent benchmark for AGN 
host comparisons, we note tha t this sample is slig htly 
larger than the one presented in Bu ndv et al.l (|2006f) be- 
cause unlike that paper, our analysis does not require 
that we restrict the perimeter of the survey on account 
of accurate environmental density measurements. While 
X-ray observations are not required for the galaxy sam- 
ple, we do select sources in the j^-band w ith good qual- 
ity ("zquality" > 3, see lDavis et al.ll2007t ) spectroscopic 
redshifts in the redshift range 0.4 < z < 1.4, reaching K s - 
band limits that depend on redshift. As in IBundv et al.l 
(j2006l), weuseK AB < 21.8, 22.0, 22.2 for z < 0.7,1.0,1.4. 

We apply the same limits to the AGN host galaxies 
which are identified in the Chandra full band and have a 
Poisson false detection probability that is less than 4 x 
10~ 6 . We then use the redshift of the matched host and 
assume T — 1.9 to infer L2-10, the X-ray luminosity 
in the 2-10 keV energy band. This effectively corrects 
for absorp tion for column densities of Nh < 10 23 cm" 2 
at z » 1 (|Nandra fe Pounds! 1 1994h . and we apply this 
procedure even if the source is not detected in the hard 
band. 

AGNs are identified as those sources with L2-10 > 10 42 
ergs s _1 ; we note that AGNs are certainly present in 
sources below this threshold, although they are more 
difficult to distinguish with the present data. Because 
normal galaxies without AGNs have not been obse rved 
with L 2 -io > 10 42 ergs s -1 fsee iBauer et al.l 12004ft , we 
adopt this limit to ensure that there is little to no con- 
tamination from sources without AGNs. Indeed, even 
without an X-ray luminosity threshold, the flux limit of 
the Chandra data (8.2 x 10 -16 ergs s _1 cm -2 in the hard 
band) suggests a very low non-AGN contamination rate 
of only a few per cent in the observed source densities 
(jBauer et al.ll2004D . 

While hard X-ray selection provides one o f the clean- 
est ways to select AGNs at high redshift (e.g.. iMushotzkvl 



12004 iBarger et al.ll2005D . it will still miss Compton-thick 
sources (Nh > 10 24 cm" 2 ), which, while not expected to 
outnumber hard X-ray detected AGNs, are thought to 
contr i bute at least half as many (e.g 



2006J iGuainazzi et"ai1 [20051 : IGilli et al.1 



Treister fc Urrvl 
20071) . We will 



account for missed Compton-thick AGNs by incorporat- 
ing the detection efficiency of our sample in the analysis 
in igl 

In what follows, we will consider the AGN host sam- 
ple with spectroscopic redshifts and Rab < 24.1 (the 
spec- z sample, 84 sources) as well as a deeper sample 
(Rab < 25.1) of AGN hosts supplemented with photo- 
m etric redshifts (the photo- z sample, 241 sources). As 
in IBundv et~aT1 (|2OO60 . photometric redshifts in AEGIS 
were estimated using two methods. For galaxies w ith 
Rab < 24.1 we use ANNz (|Collister fc Lahavl 120041 ). a 
neural network redshift estimator, which benefits from a 
large training set provided by the magnitude-limited na- 
ture of the DEIMOS target selection in the EGS. We 
also include host galaxies with 24.1 < Rab < 25.1 
and 3<r detections in t he BIK band s, but use the BPZ 
SED-based estimator (jBemtezI 120001) for these sources. 
By comparing these photo- z 's to the DEEP2 spectro- 
scopic redshifts we find Sz/(1 + z) = 0.11. Further de- 
tails on these two methods as applied to g alaxies in the 
DEEP 2 /Palomar sample can be found in IBundv et al.l 
(2006). Photometric redshifts in the EGS based on the 
CFHT Legacy Survey are also available from lllbert et al.l 
(2006). Comparing these estimates to the DEEP2 spec- 
troscopic redshifts we find 5z/(l + z) = 0.17 and there- 
fore use the slightly better ANNz+BPZ estim ates de- 
scribed above. Note that lllbert et al.l (|2006l ) redefine 



Sz/(1 + z) as 1.48 x median \z. 



spec 



•Zphotl /(I + 



'spec ) 



and that with this definition the comparison to DEEP2 
yields (fenbe r t = 0- 03, similar to the accuracy reported in 
lllbert et al.l (12001 . 

An obvious concern with photometric redshift esti- 
mates for AGN hosts is that non-thermal contamination 
could lead to large redshift errors. By comparing photo- 
metric and spectroscopic redshift estimates in the AGN 
spec-0 host sample in Figure [T] we show that in general 
this is not the case. Since the global sample of galax- 
ies (0.4 < Zphot < 1-4) is characterized by the same 
5z/(l + z) = 0.11, this comparison demonstrates that 
photometric redshifts of X-ray selected AGN hosts can be 
believed at the same level as those of non-AGN galaxies. 
This implies that AGN contaminati on in the optical/IR 
is not significant, as expected (e.g.. IBarger et al1l2005h . 
Figure [1] does show, however, that the photo- z outliers 
are dominated by the more X-ray luminous systems, sug- 
gesting more contamination in these cases. Caution must 
be used in interpreting the properties of these sources. 

4. PHYSICAL PROPERTIES 

4.1. Restframe U — B Color 

Foll owing IBundv et al.l (|2006l ) and iNandra et al 
2007), we use the methods described in IWillmer et al 



(20 
(20 



2006) to estimate the restframe (U — B) colors in both 
our galaxy and AGN host samples. This measurement is 
frequently used as a diagnostic of star formation activity 
in galax ies and exhibits a bimodal distribution to at least 
z ~ 1 (jBell et al.ll2004h . separating star- forming "blue 
cloud" galaxies from the most ly passive "red sequence." 
We used the cut employed by Ivan Dokkum et all (|2000l ) 
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Fig. 1. — Comparison of photo-z quality for the AGN spec-z host 
sample. Four additional outliers with z spcc > 2.0 are not shown, 
and all have L2— 10 > 10 43 ergs s — 1 , with three having L2— 10 > 
10 44 . As indicated, for 0.4 < z phot < 1.4, 5z/(l + z) = 0.11, the 
same as the full galaxy population. We note, however, that photo-z 
outliers tend to be the most X-ray luminous. 
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Fig. 2. — Reduced x 2 of the best fitting SED determined by the 
M, estimator as a function of restframe (U — B) color. More X-ray 
luminous AGNs tend to exhibit poorer fits, indicating the presence 
of AGN contamination. The blue colors of the AGN hosts of X- 
ray luminous sources may be caused in part by this contamination. 
The grey dashed line is typical of the division used to divide red 
and blue galaxies. 

which can be expressed in Vega magnitudes as, 



U -B = -0.032(M B + 21.52) + 0.454 - 0.25 (1) 

The (U — B) color distribution for the AGN spec-z 
host sample is illustrated in Figure [5] and was origi- 
nally discussed for a subset of the data utilized here in 
iNandra et al.1 (|2007l ) . The properties of the current and 
somewhat la rger sample here sho w agreement with those 
presented in INandra et al.l (|2007f ) . 

4.2. Stellar Mass Estimates 



To estimate s tellar masses (M*) we use the meth- 
ods described in lBundv et al.1 ([2006) . Based on the ob- 
served BRIK colors (measured using 2'.'0 diameter aper- 
tures) and the redshift information for each galaxy, we 
fit the observed SED to a grid of 13440 m odels con- 
structed using the IBruzual fe CharlotJ (|2003h population 
synthesis code. The grid spans a range of metallici- 
ties, star formation histories (parametrized as exponen- 
tial), ages, and dust content. The grid is restricted 
such that only models with ages (roughly) less than the 
cosmic age at a galaxy's redshift are considered. Sys- 
tematic uncertain ties are introdu ced by our choice of 
a Chabrier IMF (jChabrierl 120031 ) and the use of the 
iBruzual fc CharlotJ (j2T)Mr code. While some studies 
suppo rt M* estimat es based on the IBruzual fc Chariot! 
(2003) software (see iKannappan fc GawT scr 2007 ), oth- 
ers h ave pointed out potential p roblems (jMaraston et al.1 
I2006T) . In lConselice et al.l (|2007l ). we show that a prelimi- 
nary analysis with the latest Bruzual and Chariot models 
(2007, in prep.) yield mass estimates that are ^0.07 dex 
smaller on average. This issue will be discussed further 
in future work. 

With the grid defined in this way, the SED of a given 
galaxy is then compared to the model at each grid point, 
where the specific JT s -band M*/Lk ratios and val- 
ues for each model are also stored. The probability that 
each model fits the data is then summed or marginal- 
ized over the grid to yield the stellar mass probability 
distribution. The median of this distribution is taken as 
the estimate of M„, and the width provides an estimate 
of the uncertainty, typically 0.1-0.2 dex. This is added 
in quadrature to the X s -band magnitude uncertainty to 
determine the final error on M*. Stellar mass estimates 
for galaxies with only photometric redshifts suffer from 
the uncertainty in luminosity distance introduced by the 
photo-z uncertainties and the possibility of catastrophi- 
cally wrong redshift information. A study of the effect 
of photometric redshifts on stellar mass was performed 
in lBundv etall (|2005l ). 

For the case of the AGN host samples, it is impor- 
tant to consider the effects of non-thermal contamination 
which may affect the inferred stellar masses. Including an 
AGN component in the model SEDs used to estimate the 
stellar mass is not practical. Instead we investigate the 
potential error by plotting the reduced \ 2 values of the 
best fitting SED from the stellar mass grid for each mem- 
ber of the AGN spec-z host sample in Figure [2j Even 
under the best conditions, we do not require a perfect 
fit to the observed SED because our Bayesian mass es- 
timator considers a range of models when assigning the 
final mass estimate. Still, mass estimates with \ 2 > 30 
should be considered with caution. Fortunately, while 
the average \ 2 value of AGN hosts is higher than for 
galaxies without AGNs (e.g., the fraction of AGN hosts 
with x 2 > 10 is 10%, while this number for all galaxies 
is 1%) and seems to correlate with L2-107 on ly 4 hosts 
have x 2 > 30, indicating that the mass estimates for the 
AGN host sample are robust. 

5. STELLAR MASS FUNCTIONS 

5.1. Methods 

With the Af* estimates described above, we are now 
in a position to construct stellar mass functions (MFs) 
from our sample. We will describe our formalism in this 
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section and present the AGN host mass function in H5.'2{ 
The same methods are also applied to construct MFs 
for the full galaxy sample. These allow us to study the 
growing fraction of red galaxies which we will use to infer 
the star formation quenching rate in £16.11 

Constructing galaxy stellar MFs requires understand- 
ing and correcting for the limitations and compl eteness of 
the su rvey data. We adopt t he V max formalism rtSchmidtl 
119681 ) to this end, following IFundv et all (|2006h . In the 
case of the galaxy sample, the maximum volume can be 
limited by either the iT s -band depth or the Rab < 24.1 
limit used to define the DEEP2 spectroscopic sample. In 
this case, for each galaxy i in the redshift interval j, the 
value of Vmax i s given by the minimum redshift at which 
the galaxy would leave the sample, 

where dflj is the solid angle subtended by the sample 

defined by the limiting i-Q-band magnitude, K 3 Um (which 
changes depending on the redshift interval j), and dV/dz 
is the comoving volume element. The redshift limits are 
given as, 

Zhigh = TCLUl{z J max , Z 3 Kum , ZR Hm ) (3) 
Zlow = z mi n (4) 

where the redshift interval, j, is defined by [i„,z4 al ], 
z J K[ . refers to the redshift at which the galaxy would 
still be detected below the if s -band limit for that partic- 
ular redshift interval, and ZR Um is the redshift at which 
the galaxy would no longer satisfy the i?-band limit of 
Rab < 24.1. We use the best- fit SED template as deter- 
mined by the stellar mass estimator to calculate z° K , 
and ZR lim , thereby accounting for the ^-corrections nec- 
essary to compute accurate V max values (no evolutionary 
correction is applied). 

In the case of the AGN host samples, the proce- 
dure must be modified to account for the limiting X- 
ray depth. This is more complicated because the limit 
varies smoothly as a function of position within a given 
Chandra pointing. Deeper sensitivity limits correspond 
to smaller effective areas. We therefore compute Chan- 
dra sensitivity curves in the full band (corresponding to 
the selection band), accounting for the overlap with the 
Palomar near-IR data which is not complete over the 
EGS field at all depths. For a given source with X-ray 
luminosity L l x , we use a T = 1.9 power-law to estimate 
the observed flux this source would have as a function of 
redshift. We then use the sensitivity curve to compute 
the corresponding solid angle over which such a source 
could be detected as a function of redshift, d£l(z). Thus, 
we derive a second V max estimate for AGN hosts based 
on the X-ray limits: 

l-Zhigh ffl 

^max,X-ray=y dQ(z)—dz (5) 

where zi ow and Zhigh are given by the boundaries of the 
redshift interval. The final V ma x for the AGN hosts is 
taken as the smaller of the V ma x computed based on the 



R and if s -band limits and V max , x-ray ■ Typically the X- 
ray volume provides the limiting V max . 

The galaxy sample and AGN spec-z host sample make 
use of DEEP2 spectroscopic redshifts only. For this rea- 
son, additional weights must be applied to account for 
the redshift targeting selection function and success rate 
of these samples. Here w e follow the techni que described 
by IWillmer et al.l (|2006h and modified in iBundv et al.l 
((20061). Specifically, we compare the number of sources 
with good quality redshifts (zquality > 3) in a given bin 
of (B — R)/(R — I)/ Rab / -Ks-band parameter space to 
the total number of sources targ eted in that same bin . 
We adopt the "optimal" model of IWillmer et al.l (|2006h . 
which accounts for the different ways that red and blue 
galaxies are likely to be excluded from the spectroscopic 
sample. 

The situation is more complicated for AGN hosts be- 
cause applying this weighting scheme to AGN hosts as- 
sumes that the photometric sources in the correspond- 
ing color /magnitude bins also host AGNs. We there- 
fore modify the weighting scheme when it is applied to 
AGN hosts as follows. In each redshift interval we de- 
termine the ratio between the number of spectroscopic 
AGN hosts and the number of potential hosts. Poten- 
tial hosts include galaxies without X-ray detections that 
have spectroscopic redshifts, stellar masses greater than 
the completeness limit, and U — B > —0.1. This color re- 
quirement is motivated by Figure [H which demonstrates 
that most AGN hosts have such colors. For the three 
redshift intervals 0.4 < z < 0.7, 0.75 < z < 1.0, and 
1.0 < z < 1.4, we find AGN fractions of 0.08, 0.09, 
and 0.14. The AGN spec-z host weights are then de- 
termined by multiplying the "optimal" weights discussed 
above by these numbers and ensuring the weight does not 
drop below 1.0. No weighting is required for the AGN 
photo-z host samples; comparisons between the spec-z 
and photo-z AGN host samples thus provide a useful 
measure of the success of our weighting scheme. 

5.2. The AGN Host Stellar Mass Function 

We plot the AGN spec-z and photo-z supplemented 
host mass functions in Figure [3] in three redshift inter- 
vals. The green shading represents the uncertainty in the 
spec-z sample, arising primarily from number statistics. 
The corresponding mass functions of the AGN photo-z 
host sample are indicated by asterisk symbols connected 
by dotted lines. For reference, each panel also indicates 
the total spec-z AEGIS mass function (light gray shad- 
ing with solid circles), the total photo-z supplemented 
MF (triangle symbol s), and in all panels the best fitting 
MF at z ps 0.5 from IBundv et alj (|2006l ). The number 
of AGN spec-z hosts diminishes significantly in our high- 
est redshift bin. Interpretations at these redshifts rely 
on the photo-z sample only. Note that the total photo-z 
sample covers a larger area (by ~25%) than the spec-z 
sample. Thus, slight differences in the total mass func- 
tions can arise from cosmic variance, especially in the 
lowest redshift interval. 

At all redshifts, Figure [3] shows that the AGN host 
MF is roughly flat across the stellar mass range sampled. 
Comparing the AGN photo-z MFs, there is evidence that 
from z ~ 1.2 to z « 0.5 the abundance of AGN hosts 
decreases roughly by a factor of 2 at all stellar masses 
probed. Because the corresponding number density of all 
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Fig. 3. — The stellar mass function of AGN host galaxies in three redshift intervals as compared to the total galaxy stellar mass function. 
Green shading traces the MF and uncertainty of AGN hosts with spectroscopic redshifts and L2-10 > 10 42 ergs s — 1 . The asterisk symbols 
with error bars and connected by the dotted line show the AGN host MF for the photo- z supplemented sample. Total mass functions from 
the AEGIS field are shown with solid circles and grey shading (spec-z sample) and triangles (photo- z supplemented sample). The solid line 
is taken from the best fit of Bundy ct al. (2006) to the total MF at z ~ 0.5. A horizontal dashed line has been drawn at <f> = 10~ 4 in all 
panels to guide the eye. 



galaxies is lower at z > 1, the fraction of systems hosting 
AGNs increases at these epochs. This is shown explicitly 
in Figured] which plots the AGN fraction as a function of 
stellar mass. On one hand, the relatively flat MFs in Fig- 
ure [3] suggest that AGN evolution — for ex ample the de- 
clining hard X-ray luminos ity density (e.g.. IBarger et all 
120051 lHasinger et al.ll2005D — is independent of host M*. 
However, the fraction of AGN hosts shown in Figure @] 
presents a different interpretation. As galaxies continue 
assembling and their abundance grows with time, X-ray 
AGNs may be increasingly turning off — especially at the 
highest masses. This would lead to stronger evolution 
in the AGN fraction (Figure [4J accompanied by milder 
evolution in the absolute numbers of AGN hosts (Figure 

E}. 

Figure [4] also shows the result of splitting the sam- 
ple using two different X-ray luminosity thresholds. The 
dark green shading and asterisk symbols denote the AGN 
host sample with £2-10 > 10 42 ergs s _1 as in Figure 
[3l The light green shading and diamonds in Figure 0] 
provide a comparison to the host MF corresponding to 
higher X-ray luminosities of L2-10 > 10 43 ergs s" 1 . Note 
that because of th e steep decline in t he X-ray luminos- 
ity function (e.g.. iBarger et al.ll2005| ). the AEGIS sur- 
vey area is too small to effectively sample sources with 
£2-10 > 10 44 ergs s- 1 . 

While it appears that the more X-ray luminous AGNs 
are generally less abundant, we find little significant dif- 
ference in the shape of the host MF as a function of 
X-ray luminosity. Those with L 2 -io > 10 43 ergs s _1 



account for roughly one-third of the full AGN sample 
with £2-10 > 10 42 ergs s _1 , but are associated with host 
galaxies with a similar mass distribution. 

Studies of the AGN X-ray luminosity function (LF) 
show that the more luminous sources are more abun- 
dant in the past relative to the less luminous ones, 
a phenomenon often termed "AGN do wnsizing" (e.g., 
lHasinger et"aH 120051; IBarger et all 120051) . Because this 
trend is most apparent for systems brighter than the knee 
in the X-ray LF— that is AGNs with L 2 -io > 10 44 ergs 
s _1 — we would not expect the effect to be strong in this 
survey, which is too small to accurately sample such lu- 
minous AGNs. There is some suggestion, however, for 
this effect in the highest mass bin of the 0.7 < z < 1.0 
redshift interval in Figure 01 where it appears the most 
massive hosts become dominated by the brightest X-ray 
AGNs. 



6. LINKING AGNS AND QUENCHING 

We now move to the primary goal of this paper. Our 
aim is to evaluate the role of AGNs in the evolving star 
formation properties of the full galaxy population. To ac- 
complish this, we will compare the rate at which AGNs 
are triggered in galaxies of a given mass with the rate at 
which star formation is quenched at these masses. We 
will use the AGN host mass function presented in the 
previous section (coupled with the AGN lifetime) to in- 
fer the AGN trigger rate. First, however, we must char- 
acterize the mass-dependent quenching rate in the total 
population. We use the full Palomar/DEEP2 sample to 
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Fig. 4. — Fractional contribution of AGN hosts to the total 
MF in log units, shown in three redshift intervals. As in Figure 
[3] dark green shading denotes the AGN spec-z host sample with 
L2-10 > 10 42 ergs s~ x , while asterisk symbols connected by dotted 
lines denote the AGN photo-z host sample with the same L2-10 
threshold. Light green shading corresponds to AGN spec-z hosts 
with the brighter X-ray cut of L2-10 > 10 43 ergs s _1 , while dia- 
mond symbols connected by dashed lines show the corresponding 
AGN photo-z MFs. A dotted horizontal line at log / = —1.5 has 
been drawn in each panel to guide the eye. 

provide this measurement below. 

6.1. The Star Formation Quenching Rate 

We will define the quenching rate, Q, as the frac- 
tion of all galaxies in a stellar mass bin that shift to 
the red sequence per Gyr. This quantity can be de- 
rived using the methods of ^5] to plot the increasing 
fraction of red galaxies (relative to the abundance of all 
galaxies) as a function of time for various bins of stellar 
mass (Figure [5]). We will use the slope of the increas- 
ing red fraction to estimate Q. We have chosen to study 
the red galaxy fraction — as opposed to absolute num- 
ber densities — because this helps mitigate uncertainties 
caused by cosmic variance, which to first order affect the 
total n umber density meas ured in a given redshift inter- 
val fsee lBundv et ai.ll2006D . 

The use of fractional abundances also provides a better 
handle on the rate at which galaxies become red, that is 
Q. In the absence of processes that shift galaxies into 
different mass bins, quenching only alters the fraction of 
red galaxies in a given M* bin. The transfer of galax- 
ies across mass bins is constrained to be small by the 
lack of significant evolution in the shape of the total MF 
from z ~ 1. However, it is possible that both star forma- 
tion and merging may move galaxies between mass bins. 
As for merging, we make the assumption that the effect 
on the red fraction is small if merging is independent of 
galaxy type and the merging rate does not vary across 
the 0.3 dex mass bins used here. As for star formation, 
because lower mass galaxies exhibit higher SF rates, their 
evolution would tend to drive the red fraction down as 
low mass blue galaxies enter a given mass bin. However, 
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Fig. 5. — Evolving fraction of red galaxies as a function of time 
in various mass bins. The growing fractions have been fit by the 
dotted lines, excluding data points where the fraction is either 1.0 
or the data are incomplete. Excluded data points are indicated by 
the cross symbols. 



as the number of galaxies forming stars at a rate suffi- 
cient to double thei r mass over a few G yr is small for 
M* > 10 10 M Q fe.g.. lFeulner et al.ll2005l ). such an effect 
would have a small impact on the red fraction. Still, we 
emphasize that the evolving red fraction in specific mass 
bins is only an estimate of the true quenching rate. 

The buildup in the fraction of red systems is clear 
in Figure [5] and allows us to crudely fit lines to the 
fractional growth rate. In this fit, we exclude points 
at early times and low masses that fall below our ex- 
pected completeness limit as well as those points for 
which the red fraction is equal to 1.0. We take the slope 
of these lines as our estimate of Q. In the mass bins cen- 
tered at logM*/M =10.5, 10.8, 11.1, 11.4, 11.7 we find 
fractional quenching rates, Q(M»), of 8%±3%, 9%±3%, 
11%±3%, 16%±6%, and 28%±23% per Gyr. While Fig- 
ure O shows that our linear approximation adequately 
fits the data, we cannot further constrain the quenching 
rates as a function of time (or redshift). Our results are 
obviously only valid until the red fraction reaches 1.0 and 
all systems are quenched and, extrapolating the fits, we 
find that "total quenching" in these mass intervals oc- 
curs when z — -0.4, 0.0, 0.4, 0.6, 0.7, or in the assumed 
cosmology, when the cosmic age equals r age — 14.7, 11.4, 
9.3, 7.5, 7.3 Gyr. Figure[5]not only reinforces the notion 
that more massive galaxies become quenched first, but 
demonstrates the new result that they may also become 
quenched faster than their lower mass counterparts. 
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6.2. Comparing AGN Triggering and Star Formation 
Quenching 

With the star formation quenching rate measured 
above, we now utilize the AGN host mass function to de- 
rive the AGN "trigger rate", x(M*), defined as the frac- 
tion of all galaxies per Gyr in which X-ray detected AGNs 
turn on as a function of stellar mass. The rate of AGN 
triggering multiplied by the timescale over which AGNs 
are visible at X-ray wavelengths is equal to the observed 
AGN fraction (shown in Figured]). If we account for the 
AGN detection efficiency, e, of our X-ray observations, 
we can write this relation as /agn(-M*) = ex(-^*) r AGN, 
where tagn is the X-ray AGN timescale. 

This timescale is the largest uncertainty in the calcu- 
lation and must be derived from theoretical arguments. 
We will consider three estimates taken from the litera- 
ture for the average va lue of tagn- F r om the detailed 
simulations studied in I Hopkins et al.l (|2005bf ) we find 
tagn ~ 0.6 Gyr. From the statisti cal and populat i on ar - 
guments for low efficiency AGNs in lMarconi et al.l (2004) 
we use tagn ~ - 9 Gyr , and from the model discussed 
in iGranato et all (|2004h we use tagn » 1.8 Gyr. We 
will return to the timescale problem and discuss these 
estimates further below. 

Assuming e = 1 we solve for the corresponding trig- 
ger rates, x(M*), by dividing the AGN host fraction 
(L2-10 > 10 , Figured]) by the timescales above. We 
plot the results for the two redshift bins where our sam- 
ple is most complete in Figure |6] The corresponding 
star formation quenching rates from our analysis of red 
galaxies is denoted by the red shaded region. The ef- 
fect of a lower detection efficiency of e w 0.7 (as might 
arise from missed Compton-thick sources) increases the 
trigger rates as roughly shown by the arrow. 

Despite the uncertainties and assumptions, Figure [6] 
demonstrates surprising agreement in both the normal- 
ization and mass dependence of the rates of quenching 
and AGN triggering, given the three estimates for tagn- 
We interpret this as strong but circumstantial evidence 
that the quenching of star formation and AGN activity 
are physically related. We will turn to the question of 
whether AGNs actually cause quenching in £16.41 

6.3. The X-ray AGN Timescale 

Because the calculation of the AGN trigger rates shown 
in Figure[6]relies heavily on the assumed value of tagnj in 
this section we further explore the X-ray AGN timescale 
and the reliability of the estimates we have used. As a 
point of reference, we begin by deriving the value of tagn 
that would be necessary to force the observed quenching 
and triggering rates to be equal. We set Q(M*) = x(M») 
and solve for tagn at each stellar mass bin where esti- 
mates of the two rates are available (essentially we divide 
the full AGN host fraction in Figure [4] by the quenching 
rates derived in Figured]). In principle tagn may be re- 
lated to the host stellar mass, but we will ignore this and 
average over the values of tagn derived for each mass 
bin to roughly estimate the range or "probability dis- 
tribution" of timescales needed to perfectly match the 
quenching and AGN trigger rates observed in our two 
redshift intervals. These are plotted as the solid and 
dashed red lines in Figure [7] and illustrate where theo- 
retical estimates of tagn would fall if it were true that 
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Fig. 6. — Comparison in two redshift bins between the frac- 
tional AGN "trigger rates" — calculated using three estimates of 
the X-ray AGN timescale — and the star formation quenching rate 
denoted by the red shaded region. The width of the shading 
illustrates the lcr uncertainty. Trigger rates are derived based 
on the AGN photo-z sample (similar results are obtained for 
the spec-z sample), assuming all AGNs are detected and X-ray 
AGN timesc ales of 0.6, 0.9, and 1.8 Gyr based on es t imates from 
the work of IHopkins et al.l H2005bl) , IMarconi efakl J200l) , and 
IGranato et al.l H2004T) . The systematic effect of a 70% detection 
efficiency (e as 0.7) is shown by the arrow 



Q(M„) = x(A'Q. 

We can now compare the estimates of tagn we have 
used and discuss their uncertainties. The most appro- 
priate predictions come from detailed s imulations an- 
alyzed by Hopkins and collaborators. Hopki ns et al.1 
(2005b) conduct five hydrodynamical simulations of gas 
rich mergers of disk galaxies that host super massive 
black holes (SMBHs). During the simulations, gas be- 
comes funneled to the center of the system, fueling the 
growth of the newly merged SMBH. The authors use 
a prescription in which some fraction (e r = 0.1) of this 
accreted material is radiated in a quasar phase. They as- 
sume 5% of this energy couples to the surrounding gas, 
helping to regulate the flow. The X-ray luminosity (and 
column density) as a function of time is calculated by as- 
suming a quasar continuum SED and ray-tracing many 
lines of sight through the gas and dust in the simulated 
galaxy. Their Figure 2 presents a relation between the 
observed AGN lifetime (tagn) and its X-ray luminos- 
ity. By applying this relation to the X-ray luminosities 
observed in our sample we derive a rough distribution 
of predicted AGN timescales which we plot as the dot- 
ted line in Figure fT] Note that the predicted tagn from 
the IHopkins et al.1 (|2005b[ ) models increases rapidly for 
AGNs with Lx < 5x 10 42 ergs s _1 , accounting for the tail 
towards longer lifetimes shown in the figure and yielding 
the estimate of tagn ~ 0.6 Gyr used in Figure [U 
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Fig. 7. — Range of X-ray detectable AGN timescales. As a 
reference point, the result of assuming the AGN trigger rate equals 
the quenching rate yields timescales in the range indicated by the 
red solid (z ~ 0.5) and dashed (z ~ 0.9) distributions. The effect 
of a 70% AGN detection efficiency (e = 0.7) on these distributions 



is shown by the re d arrow. Independent p redictions of tagn based 
on the models of Hopkins et al. (2005b) have a range indicated 
by the dotted line. Predictions from IGranato et alj II2004D and 
Marconi et al. (2004) arc denoted by the grey shaded regions. 



The work of IGranato et alj (|2004f ) provides another 
independent comparison, also based on numerical sim- 
ulations that encode the effects of star formation, cool- 
ing, supernovae feedback, and AGN feedback set in the 
co ntext of dark matter halos. The systems analyzed 
in IGranato et al.l (|2004f ) are non-interacting spheroidal 
galaxies as opposed to merging disks, however, and while 
the simulations lack detailed modelling of the AGN X- 
ray emission, some rough constraints can be obtained for 
the predicted values of tagn- Their Figure 3 shows the 
black hole accretion rate as a function of cosmic time in 
their simulations. Since, as we show below, our AGNs 
are likely accreting at significa ntly sub- Eddington r ates, 
a rough estimate of tagn from lGranato et all (|2004T ) can 
be gained by measuring the time between peak black hole 
accretion in their simulations and the point at which the 
accretion drops below a factor 0.01-0.001 of the maxi- 
mum rate. This yields AGN timescales between ~1 and 
~2 Gyr as shown in Figur e [71 by the grey shade d region. 
For an average value from IGranato et al.l (|2004f ) we take 
Tagn ~ 1-8 Gyr. 

Finally, much work has been invested in utilizing vari- 
ous observations to constrain the lifetime of bright quasar 
activity, typic ally resulting in values of 10 7 -10 8 yr (see 
the review by iMartrnl |2004[) . However, historically the 
focus has rested on the brightest quasar phase during 
which black hole growth is thought to be most rapid. 
This phase does not correspond to the l ower accretion 
rates o f the AGNs detected in our sample. iMarconi et al.l 
(|2004h . however, provide a suitable estimate based on 
matching the local black hole density to the AGN lumi- 
nosity function. For low efficiency sources, the predicted 

range is tagn ~ 0.5 1 Gyr, weighted more towards 

1 Gyr for the lowest luminosity AGNs and providing a 
rough average value of tagn = 0.9 Gyr. 

It is clear from Figure [JJ that la rge uncertainties exist . 
Still, from the detailed analysis in lHopkins et al.l (l2005bf) 



to the mor e approximate estima tes from IMarconi et alJ 
( 2004) and IGranato et all (|2004h . this plot demonstrates 



that the predicted range and uncertainty in the timescale 
for X-ray AGN activity using a variety of methods is at 
least compatible with a scenario in which AGN trigger- 
ing is linked to quenching. This helps to validate the 
agreement seen in Figure [6l Clearly further progress in 
confirming this link would strongly benefit from addi- 
tional detailed predictions of the X-ray properties and 
lifetimes of AGNs. 

6.4. Is AGN Feedback Responsible For Quenching? 

We have argued that the similarity in the rate of AGN 
triggering compared to the rate of star formation quench- 
ing suggests that the two phenomena are linked. But 
what is the nature of this link? Specifically, we would 
like to know if AGNs, perhaps through feedback mecha- 
nisms, are directly responsible for the quenching of star 
formation. In this section we will begin to probe this 
question by investigating the individual properties of our 
X-ray selected AGNs and how they correlate with their 
host galaxies. 

Our strategy will be to study the SMBH accretion 
rates — as parametrized by the Eddington ratio — in our 
sample. We acknowledge that these estimates are some- 
what crude and subject to systematics, but argue that 
they nonetheless provide important insight on how AGNs 
are related to the properties of their host galaxies. The 
Eddington ratio compares the bolometric luminosity, 
Lbolj of the AGN to the Eddington Luminosity which 
is simply related to the SMBH mass, ^Edd = 1-25 x 
1O 38 (M BH /M ) ergs s _1 . To determine M B h we use the 
relation between if- band bulg e luminosity and black hole 
mass calcula ted bv iGrahaml (I2007D and based on pre- 
vious work (Marconi k Huntjl^ol iMcLure fc Dunlopl 
I2004h . Note that IWoo et alJ (|20Q6h find evidence that 
the z ~ 0.3 relation is offset relative to that at z — by 
a +0.6 dex increase in Mbh, although we do not apply 
this correction here. As most X-ray selected AGN hosts 
at th e redshifts of our sample have early- typ e morpholo- 
gies (jGrogin et al.l 120051 : iPierce et al.l l2007p . we assume 
that the bulge-to-total ratio is 1.0, but in what follows we 
will demonstrate the effect of lower ratios on our results. 
Finally, we calculate I/bo i by assuming a h ar d X-ray bolo- 
metric correction of 35 (jElvis et al.lll994D . iBarger et al.1 
(2005) argue that for obscured ( narrow line) AGNs , the 
correction should be 85 although iPozzi et"aTI (|2007l ). us- 
ing a sample of type-2 AGNs, find a wide range of bolo- 
metric corrections with an average of ^25. Clearly the 
bolometric corrections are uncertain within a factor of 
-2-3. 

Using the methods above and with the stated caveats 
in mind we plot the estimated Eddington ratio as a func- 
tion of stellar mass in Figure [8l The X-ray luminosity 
range of the data is indicated, showing that more lumi- 
nous AGNs tend to have more efficient accretion. The 
isolated error bar indicates the typical uncertainties from 
observational scatter, while the arrows show the effects of 
systematic errors in estimating the bulge luminosity and 
bolometric correction. The "Bulge" systematic shows 
what would happen if the bulge-to-total ratios were low- 
ered by a factor of 2, or equivalently if the applied value 
of AfeH/A^buige was increased by the same amount. The 
"Bolometric" systematic shows the effect of increasing 
the bolometric correction from 35 to 85. It is important 
to note that our calculation of the Eddington ratio en- 
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Fig. 8. — Eddington ratios in the AEGIS spec-z sample plotted 
against the host galaxy stellar mass. Different symbols indicate dif- 
ferent ranges of AGN X-ray luminosity as indicated. The isolated 
error bar illustrates the typical uncertainty arising from observa- 
tional scatter. The labelled arrows suggest the magnitude and 
sense of systematic uncertainties. The "Bulge" systematic shows 
what would happen if the bulge-to-total ratios were lowered by a 
factor of 2, or equivalently if the applied value of iWBjj/A^bulge 
was increased by the same amount. The "Bolometric" system- 
atic shows the effect of increasing the bolometric correction from 
35 to 85. Note that the axes are not fully independent because 
the Eddington ratio is proportional to M^T , leading to the slight 
downward trend observed. 

sures that it is proportional to M ~ 1 so that downward 
mass-dependent trends in Figure [5] are expected given 
our methodology. 

What is perhaps most revealing about the figure, how- 
ever, is the large range (more than 2 orders of magnitude) 
in Eddington ratios that is apparent for host galaxies of 
all masses. Figure [3] shows a similar diagram where it 
is now possible to compare the Eddington ratios versus 
host galaxy restframe (U — B) color. The red and blue 
distributions of all galaxies are indicated by the back- 
ground shading. Note the existence of very blue hosts 
(beyond the range of normal c olors) dominated by b right 
X-ray systems. As observed in lNandra et al.l (|2007f ). it is 
likely that the AGN contaminates the host color of these 
galaxies. 

Considering the hosts in the blue cloud and red se- 
quence, there appears to be little difference in the large 
spread of accretion rates. This is especially the case when 
one ignores the X-ray brightest systems (marked as tri- 
angles), whose true host colors may be redder than ob- 
served. We note that among all galaxies at z « the 
fraction of dusty but still star-forming systems on the 
red sequence is ~ 7% and is not li kely to be significantly 
higher at z ~ 1 ()Yan et al.ll2006f h It is therefore likely 
that most red AGN hosts in Figure [Dare truly quenched 
systems. 

While the uncertainties in Figures [8] and [9] are large, 
the intrinsic scatter appears to be more substantial and 
suggests little or no additional trends with stellar mass 
or color beyond those expected from the methodology. 
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Fig. 9. — Eddington ratios as in Figure [8] for host galaxies of 
varying restframe (U — B) color. A typical value for the division 
between red and blue galaxies is shown with the vertical dotted 
line and the color distributions of all galaxies is indicated by the 
background blue and red shading. 

These observations may therefore have important impli- 
cations for the question of whether AGN feedback is re- 
sponsible for quenching. 

Some feedback models imagine that the AGN has an 
explosive episode that drives gas out of the galaxy halo. 
In the model of Hopkins and collaborators discussed 
above, for example, AGN activity is triggered by major 
mergers. The buried AGN is virtually undetectable until 
a violent quasar phase in which tremendous energy is ex- 
pelled, sufficient to heat or dispel most of the remaining 
gas in the gala xy, thus quench i ng furt her star formation. 
As described in lHopkins et al.l (|2007bf ). in this model one 
would expect that AGN activity evolves with time as 
AGN feedback during the quasar phase impacts the sur- 
rounding material, heating and driving it from the newly 
merged system. The X-ray observations studied here 
would then correspond to the post-quasar phase tracing 
AGNs as they decay to low luminosities; beforehand the 
obscuration is predicted to be >10 24 cm -2 , enough to 
absorb hard X-rays. The short-lived quasar phase that 
immediately follows w ould be too bright to en able studies 
of host properties (see iHopkins "et~al1l2007bl) . 

Figures [8] and [9] suggest some potential problems with 
simple interpretations of explosive models of this sort. 
One might expect the most massive and reddest hosts 
to have been quenched earliest and therefore to harbor 
AGNs in the latest stages of decay. Little fuel should 
remain in such systems long after the quasar phase, 
and the luminosities and Eddington ratios should be 
low. This picture appears consistent with observations 
of local AGN hosts identified through optical emission 
line diagnostics. Stronger emission line AGNs are found 
in younger stellar populations with higher specific star 
formation rates, while weaker AGNs favor hosts that 
have apparently b een quenched (Kauffm ann et al"1l2003h 
ISalim et aill2007t ). Figures [5] and H show that X-ray se- 
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lected AGNs may present a different picture. Here, the 
reddest and most massive hosts harbor AGNs that cover 
nearly the full range of X-ray luminosity and accretion 
rates. Indeed, Figure [3] demonstrates that AGNs are 
hosted by galaxies with masses covering the full range 
probed. 

One perspective on this quest ion is suggested by the 
work of ICiotti k, Ostrikerl (|2007l ) who demonstrate that 
X-ray luminous AGN activity as well as starbursts can be 
effectively fueled by stellar mass loss from evolved stars 
in old stellar populations. While it is not clear what 
timescales would be involved, in the absence of other 
fueling mechanisms, this process requires an old stellar 
population to function, generating an obvious link be- 
tw een quenching and the ap pearance of AGNs. The work 
of ICiotti &; Ostrikerl (|2007l ) serves to demonstrate that 
AGNs may be fueled by a variety of mechanisms. Even 
if a high-accretion quasar phase was initially responsible 
for disrupting the internal gas supply, the AGN may later 
be "refueled" by such mechanisms, including the inflow 
of gas lost from evolved stars. 

Finally, in addition to refueling, an alternative expla- 
nation for the range of AGN properties observed in the 
most massive and reddest hosts could come from the no- 
tion that AGN feedback is not responsible for quench- 
ing star formation, but is triggered by the same process 
that is. A nuclear starburst fueled by the same inflowing 
gas that ignites the AGN is a promising example, capa- 
ble of providing the feedback energy necessary (in the 
form of stellar winds) to heat and expel the surround- 
ing gas supply and help regul ate correlations betwee n 
bulge and SMBH properties (see lFerrarese fc Fordl feOOS). 
Previous work has revealed evidence for a connection 
between starbursts and AGN activit y (e.g.. | Yan et alJ 
2006t iGotol 12001 lYa^ elall [200& IWild et all l2007t 
Georgakakis et al.ll2007al) which has also been explored 
in mo dels (e.g., ISomerville et al.l 120011 : I Hopkins et alJ 
l2007ah . 

A mixed AGN/star burst scenario is supported by the 
observations presented here. AGN activity could be ini- 
tially triggered during or towards the end of starburst 
quenching, leading to a greater diversity in the phases 
of X-ray detecte d AGN activi t y amo ng galaxy hosts. 
As suggested by ICroton et alJ (|2006f h low luminosity 
AGNs — undetected in the curren t sample — may be ubiq- 
uitous in quenched systems (e.g.. lSalim et al.ll2007l ). pro- 
viding the necessary feedback that prevents further star 
formation. As part of this feedback cycle, these systems 
could periodically enter active phases that could be de- 
tected through X-ray emission. 

7. SUMMARY 

We have used the combination of the DEEP2 Galaxy 
Redshift Survey, Palomar near-IR imaging, and Chandra 
X-ray observations to study the properties of galaxies 
that host X-ray selected AGNs. We summarize our find- 
ings below. 

• The AGN host stellar mass function over the red- 
shift range 0.4 < z < 1.4 is roughly flat as a 
function of M* . The abundance of AGNs appears 
higher at z ~ 1 by a factor of ~2. Coupled with 
the decrease in number density of all galaxies at 



high redshift, the AGN fraction increases at early 
times, especially among the most massive galaxies. 

• The MF of host galaxies for an X-ray luminous 
subset of our sample with AGN luminosities of 
£2-10 > 10 43 ergs s _1 indicates a lower abundance 
but a very similar mass dependence as the full sam- 
ple with L2-10 > 10 42 ergs s _1 . We see some evi- 
dence, however, that brighter AGNs, whose abun- 
dance increases with redshift, are hosted by more 
massive galaxies. 

• Using the full DEEP2/Palomar sample we estimate 
the star formation quenching rate, defined as the 
number of galaxies that move to the red sequence 
per Gyr. Our estimates suggest that massive galax- 
ies not only populated the red sequence at earlier 
epochs but did so at a faster rate than less massive 
galaxies. 

• We show that the quenching rate agrees with the 
rate at which AGN activity is triggered in galaxies 
if the lifetime over which AGNs would be detected 
through X-ray emission is 0.5-1 Gyr, similar to es- 
timates from detailed predictions based on numer- 
ical simulations and black hole demographics. The 
agreement between the mass-dependent quenching 
and AGN triggering rates is evidence of a physical 
link between these two phenomena. 

• We test the causality of this link by comparing 
black hole accretion rates to the stellar mass and 
color of associated host galaxies. The most massive 
and red hosts — which presumably have quenched 
at the earliest times — harbor X-ray selected AGNs 
as active as those found in blue hosts. This suggests 
a more complicated relationship between AGNs 
and star formation. It is possible that X-ray se- 
lected AGNs are associated with but do not directly 
cause star formation quenching and, furthermore, 
may be subject to refueling after quenching occurs. 



We would like to thank Ray Carlberg for useful dis- 
cussions on this project. KB would like to acknowledge 
support from Olivier Le Fevre and the Observatoire As- 
tronomique de Marseille Provence where he was a visiting 
researcher. The Palomar Survey was supported by NSF 
grant AST-0307859 and NASA STScI grant HST-AR- 
09920. 01-A. Support from National Science Foundation 
grants 00-71198 to UCSC and AST 00-71048 to UCB 
is also gratefully acknowledged. Financial support has 
also been provided through PPARC and the Marie Curie 
Fellowship grant MEIF-CT-2005-025108 (AG), the Lev- 
crhulme trust (KN), the Hubble Fellowship grants HF- 
01165. 01-A (JAN) and HF-01182.01-A (ALC) and the 
STFC (EL). We wish to recognize and acknowledge the 
highly significant cultural role and reverence that the 
summit of Mauna Kea has always had within the in- 
digenous Hawaiian community. It is a privilege to be 
given the opportunity to conduct observations from this 
mountain. 



12 



Bundy et al. 



REFERENCES 



Barger, A. J., Cowie, L. L., Mushotzky, R. F., Yang, Y., Wang, 
W.-H., Steffen, A. T., & Capak, P. 2005, AJ, 129, 578 

Bauer, F. E., Alexander, D. M., Brandt, W. N., Schneider, D. P., 
Treister, E., Hornschemeier, A. E., & Garmire, G. P. 2004, AJ, 
128, 2048 

Bell, E. F., Wolf, C, Meisenheimer, K., Rix, H.-W., Borch, A., Dye, 
S., Kleinheinrich, M., Wisotzki, L., & Mcintosh, D. H. 2004, ApJ, 
608, 752 

Bem'tez, N. 2000, ApJ, 536, 571 

Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Borch, A. et al. 2006, A&A, 453, 869 

Bower, R. G. et al. 2006, MNRAS, 370, 645 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Bundy, K., Ellis, R. S., & Conselice, C. J. 2005, ApJ, 625, 621 

Bundy, K. et al. 2006, ApJ, 651, 120 

Canalizo, G. & Stockton, A. 2001, ApJ, 555, 719 

Cattaneo, A., Dekel, A., Devriendt, J., Guiderdoni, B., & Blaizot, 
J. 2006, MNRAS, 370, 1651 

Chabrier, G. 2003, PASP, 115, 763 

Cimatti, A., Daddi, E., & Renzini, A. 2006, A&A, 453, L29 

Ciotti, L. & Ostriker, J. P. 2007, ApJ, 665, 1038 

Coil, A. L., Newman, J. A., Kaiser, N., Davis, M., Ma, C.-P., 

Kocevski, D. D., & Koo, D. C. 2004, ApJ, 617, 765 
Collister, A. A. & Lahav, O. 2004, PASP, 116, 345 
Conselice, C. J. et al. 2007, preprint (arXiv0708.1040), 708 
Cooper, M. C. et al. 2007, preprint (arXiv0706.4089), 706 
Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 

112, 839 

Croton, D. J. et al. 2006, MNRAS, 365, 11 

Cuillandre, J.-C, Luppino, G., Starr, B., & Isani, S. 2001, in SF2A- 

2001: Semaine de TAstrophysique Francaise, 605 — h 
Davis, M., Gerke, B. F., Newman, J. A., & the Deep2 Team. 2005, 

in ASP Conf. Ser. 339: Observing Dark Energy, ed. S. C. Wolff 

& T. R. Lauer, 128-+ 
Davis, M. et al. 2003, in Discoveries and Research Prospects from 

6- to 10-Meter-Class Telescopes II. Edited by Guhathakurta, 

Puragra. Proceedings of the SPIE, Volume 4834, pp. 161-172 

(2003)., 161-172 
Davis, M. et al. 2007, ApJ, 660, LI 
de Lucia, G. & Blaizot, J. 2007, MNRAS, 375, 2 
Dunlop, J. S., McLure, R. J., Kukula, M. J., Baum, S. A., O'Dea, 

C. P., & Hughes, D. H. 2003, MNRAS, 340, 1095 
Elvis, M. et al. 1994, ApJS, 95, 1 
Faber, S. M. et al. 2007, ApJ, 665, 265 

Ferrarese, L. & Ford, H. 2005, Space Science Reviews, 116, 523 
Feulner, G., Gabasch, A., Salvato, M., Drory, N., Hopp, U., & 

Bender, R. 2005, ApJ, 633, L9 
Georgakakis, A. et al. 2006, MNRAS, 371, 221 
— . 2007a, in prep., 660, L15 
— . 2007b, ApJ, 660, L15 
Gerke, B. F. et al. 2007, MNRAS, 376, 1425 
Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 
Goto, T. 2006, MNRAS, 369, 1765 
Graham, A. W. 2007, MNRAS, 543 

Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 

2004, ApJ, 600, 580 
Grogin, N. A. et al. 2005, ApJ, 627, L97 

Guainazzi, M., Matt, G., & Perola, G. C. 2005, A&A, 444, 119 



Hasinger, G., Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417 
Hopkins, P. F., Cox, T. J., Keres, D., & Hernquist, L. 2007a, 

preprint (arXiv0706.1246), 706 
Hopkins, P. F., Hernquist, L., Cox, T. J., & Keres, D. 2007b, 

preprint (arXiv0706.1243), 706 
Hopkins, P. F., Hernquist, L., Martini, P., Cox, T. J., Robertson, 

B., Di Matteo, T., & Springel, V. 2005a, ApJ, 625, L71 
Hopkins, P. F. et al. 2005b, ApJ, 630, 705 

Hutchings, J. B., Cherniawsky, A., Cutri, R. M., & Nelson, B. O. 

2006, AJ, 131, 680 
Ilbert, O. et al. 2006, A&A, 457, 841 
Juneau, S. et al. 2005, ApJ, 619, L135 
Kannappan, S. J. & Gawiser, E. 2007, ApJ, 657, L5 
Kauffmann, G. et al. 2003, MNRAS, 346, 1055 
Maraston, C. et al. 2006, ApJ, 652, 85 
Marconi, A. & Hunt, L. K. 2003, ApJ, 589, L21 
Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R., & 

Salvati, M. 2004, MNRAS, 351, 169 
Martin, D. C. et al. 2007, preprint l |astro-p h/0703281 1 
Martini, P. 2004, in Coevolution of Black Holes and Galaxies, ed. 

L. C. Ho, 169-+ 
McLure, R. J. & Dunlop, J. S. 2004, MNRAS, 352, 1390 
Mushotzky, R. 2004, in Astrophysics and Space Science Library, 

Vol. 308, Astrophysics and Space Science Library, ed. A. J. 

Barger, 53 — h 

Nandra, K. & Pounds, K. A. 1994, MNRAS, 268, 405 

Nandra, K. et al. 2005, MNRAS, 356, 568 

— . 2007, ApJ, 660, Lll 

Pierce, C. M. et al. 2007, ApJ, 660, L19 

Pozzi, F. et al. 2007, A&A, 468, 603 

Salim, S. et al. 2007, preprint (0704.3611), 704 

Scannapieco, E., Silk, J., & Bouwcns, R. 2005, ApJ, 635, L13 

Schmidt, M. 1968, ApJ, 151, 393 

Silk, J. & Rees, M. J. 1998, A&A, 331, LI 

Somerville, R. S., Primack, J. R., & Faber, S. M. 2001, MNRAS, 
320, 504 

Treister, E. & Urry, C. M. 2006, ApJ, 652, L79 
Treu, T., Ellis, R. S., Liao, T. X., & van Dokkum, P. G. 2005, ApJ, 
622, L5 

van Dokkum, P. G., Franx, M., Fabricant, D., Illingworth, G. D., 

& Kelson, D. D. 2000, ApJ, 541, 95 
Wild, V., Kauffmann, G., Heckman, T., Chariot, S., Lemson, G., 

Brinchmann, J., Reichard, T., & Pasquali, A. 2007, preprint 

(arXiv0706.3113), 706 
Willmer, C. N. A. et al. 2006, ApJ, 647, 853 

Wilson, J. C. et al. 2003, in Instrument Design and Performance 
for Optical/Infrared Ground-based Telescopes. Edited by lye, 
Masanori; Moorwood, Alan F. M. Proceedings of the SPIE, 
Volume 4841, pp. 451-458 (2003)., 451-458 

Woo, J.-H., Treu, T., Malkan, M. A., & Blandford, R. D. 2006, 
ApJ, 645, 900 

Yan, R., Newman, J. A., Faber, S. M., Konidaris, N., Koo, D., & 

Davis, M. 2006, ApJ, 648, 281 
Yang, Y., Tremonti, C. A., Zabludoff, A. I., & Zaritsky, D. 2006, 

ApJ, 646, L33 



